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PREFACE 


The operation of a nuclear power plant involves a variety of 
instrumentation, most of which is of a conventional character. However, 
for control purposes, conventional instrumentation for power level 
determination is not adequate since it does not respond instantaneously 
either to overall or to local changes in reactor power. The reactor 
power level is proportional to the neutron flux, and instruments that 
measure flux can be located at many points throughout the reactor system. 
Such devices have an essentially instantaneous response and are suitable 
for indicating power levels and for operating automatic control and safety 
mechanisms, 

With respect to power level determination, the Westinghouse Corpora- 
tion, through operating experience with the prototype of the Submarine 
Thermal Reactor, has determined that there are two basic deficiencies in 
the neutron detection systems. These are that (a) an array of detection 
instruments is required in order to minimize the effect of changes in 
control rod configuration, thus complicating the instrumentation by 
requiring either multiple indicators or somes sort of scanning-averaging 
circuit, and (b) there is an aging phenomenon associated with neutron 
detection instruments when used in high flux regions for extended pericds, 
thus requiring periodic calibrations with thermal instrumentation and 


eventually replacement, 
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It would be desirable to have a system for power level determination 
which incorporates the features of the neutron detection systems, without 
these deficiencies, It is the aim of this paper to describe one such 
System, and evaluate its probable behavior under different operating 
conditions. 

In order to avoid any breach of security, the writer has not used 
any actual dimensions or characteristics which would apply specifically 
to the Submarine Thermal Reactor. However, an attempt has been made to 
arrive at results which are close enough to the actual case to have some 
meaning when compared with the experimental results. For this same reason, 
any description of actual systems is either in very general terms or left 
to the reader to look up in the classified material available. 

Experimental work in connection with this paper was done at the 
Submarine Thermal Reactor Test Site, Arco, Idaho, with the assistance 
and cooperation of the employees of the Westinghouse Corporation at that 
facility. The writer particularly wishes to thank Doctor W. H. Esselman, 
Head of the Technical Section, who presented the original idea for 
investigation, and Mr. Grady B. Matheny of that section who served as a 
technical advisor during the experimental work. Also, Professor E. E. 
Drucker of the U. S. Naval Postgraduate School, for his assistance, 


encouragement, and cooperation in the preparation of this paper. 
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CHAPTER I 


INTRODUCTION 


1. If N radioactive nuclei are present at time t, and if no new nuclei 
are introduced into the sample, the rate of change of active nuclei with 


respect to time is 


pir > -AN (1) 
; mS 
“ NM etNe (2) 
where A =the disintegration constant — 0.693 (3) 
T 
a or something proportional to git is observed rather than N. 


Differentiating (2) and substituting from (1) 
dN _sdN) _-ak 
dt era = 
or R= R, e** (4) 
The disintegration rate R is called the activity of the sample. 
2. In chain decay equation (2) holds for the first member. Let 
species A decay into species B with disintegration constant A, 3; Let 
B decay into C with disintegration constant Ag ; one does not care 
what happens to C. Let Na and Ne be the number of A and B atoms present 
at time t. For Ny one has from equation (2) 
-A,t 
= Nye (5) 
dng 
The net rate at which Np is changing (3 ) is the difference 


between the rate at which it is being produced (- so) and the rate 


N, 


of which it is decaying [A,%; See equation (1) ] . 


Using equation (5) this gives 


dA Nz 


-A,t 
ge = AaNaoe ” 


Ag N, (6) 





Solving this differential equation for N, and assuming that one 


B 
starts with N, = 0 at t-0 leads to 





= An - Aat Age 
NG hee. Wao (2 - © ) (7) 
3. A limiting case of equation (7) is that in which A is very long-lived 


-Apt 
or A,<o AQ- At large t, e7*8t can always be neglected relative to e “, 


Np then approaches 
= Awe ; * 
eo Naan © Tig) (8) 


é 


If the period of A is not only much larger than that of B but also is 
large compared with the time during which one makes observations, ep Xat 
will not change much during the observation, and N, will approach a 
constant "secular equilibrium" value N Boo given by 
is We el (9) 
This same value of NB wo is approached no matter what the initial 
value of Ne is, 
Since A,<cA, with An = O , and substituting from equation (9), 
equation (7) may now be written . 
Ne = No (lee*#*) (10) 
Secular equilibrium occurs when a radioactive substance is being 
made at a constant rate either by the decay of its long-lived parent 
(as above) or by bombardment in a cyclatron or nuclear reactor, In 
either event, after waiting approximately six half-lives of B, Ne would 
approach its equilibrium value, where the decay rate (N,., A, ) equals 
the constant production rate (N, An )g See equation (9). 
4, For the particular case under consideration (i.e. neutron reactions 
with the coolant in a reactor producing a gamma radioactive species) the 
constant production rate is xq, where £ is the macroscopic cross-section 


of the coolant and ¢ is the flux within the reactor core (which will be 


assumed uniform and monoenergetic for purposes of this paper). Sub- 





stituting this value for N, A, in equation (9) leads to the secular 


equilibrium value given by 


£¢ 
N a0 ae 
or Nu, A, = 29 (11) 


If t = 0 is considered to be the time a unit volume of coolant 
leaves the region of high flux; the instantaneous initial decay rate, 
by substitution of equation (11) in equation (1), is found to be 

c YY, = Ni ay 2 £9 (12) 
At a time t after leaving the high flux region the instantaneous 
decay rate will be, from equation (4) 
ANe _ £Ms ) on ret 
dt dt fo 
or R = R,e748* (13) 
Where R, is proportional to £ ¢@. 
(i.e. R HK'E¢ ) (14) 

It may occur to the reader that the unit volume of coolant leaves 
the reactor core, traverses the loop piping, then returns to the high 
flux region, and therefore may not actually reach the secular equilibrium 
condition before leaving the reactor again. This is not the case, as the 
relative magnitudes of the times involved (i.e. loop transit time, time 
within the core, and half life) are such as to minimize this effect. 
This will be discussed further in the next chapter using the system 
described there for examples. 

5. The power level of a given power reactor is obviously directly 
related to the rate of heat release within the core. The rate of heat 


release is dependent on the fission rate, The fission rate for the 


reactor is = q@ , therefore 
fuel 


=_ i! 
Power =K 2. ¢g (15) 





Since for a given reactor installation K" and 2 will be 
fuel 
constant, power is directly proportional to flux (or vice versa). 


Substituting équgt doris ((14) iand1@)s5) inte equation (13) gives 


= =i t ’ -X t Po w -Aae 
R= Roe? swe gie a(n’ (ES \} 
Fuel 
or R= KP aot (16) 
where K = fee 
K Ze igl 


up Percent Power 
06 
6. The time, t, in equation (16) will be the time required for a given 
unit volume to travel from the reactor core outlet to a fixed point of 
measurement. 
Let A = Flow Area of Coolant Loop Piping 
F = Fraction of Full Flow (as a decimal) 


Length of Piping from Core Outlet to 


L 


Point of Measurement 


# 
if 


Capacity of Coolant System 


3 
il 


Time for a Unit Volume to Complete 
the Flow‘Cycle at Full Flow 


then 
TAL 
t= OF | (17) 


Substituting in equation (15) 


TAL 
R = KP eel ar (18) 


for a given system this may be written as 
-S$ 
R= KP e * (19) 


Ag TAL 
Q 


where C = 





(3) N'(n Res, nl7 0 yl? emits a 3.7 Mev beta in going to 
ol? then a 0.9 Mev neutron and ends 
as or 
3. It can be seen from the preceeding section that there are only four 
possible reactions which should be considered. These are the 0(7,P) 
ni6 / ol7(n p) yi? , ol8(n x) ol 9 , and the nt®( n+ ) ni? Using 


Feathers formula for a rough range—energy relation for the beta radiation [4] 


range (g/cm*) = 0,546E(Mev) — 0.016 (1) 
one finds the range of the most energetic beta emitted, from nl to be 
range = (0.546)(10) — 0.016 = 5.30 g/cm*. (2) 


With an average density of steel, used for reactor structure and piping 
assumed as 7.9 g/cm, this leads to a thickness of steel of less than 
one cm, = = 0.67) to reduce this radiation to a point where it is 
indistinguishable from the background. It is therefore reasonable to 
assume that any beta activity can be neglected. 

The threshold for the 017 (n,p)Nn2? reaction, which results in the 
emission of the 0.9 Mev neutron, is 7.9 Mev, The average cross section 
above the threshold is reported as 0.0] barn for an average neutron 
spectrum well above the threshold energy, but is not considered more 
reliable than by about a factor of 10, The natural abundance of 01? 
is 0.039%. Considering the abundance, reaction cross-section, and 
the attenuation probability, it is reasonable to assume that any effect 
from this neutron emission will be negligéble. 

The isotope o18 has an abundance of 0.20% and a thermal capture 
cross-section of 2.2 x 10" barn, The 018 (n, ¥ )ol9 reaction emits a 


1.6 Mev delayed gamma with a 29 second half-life. Comparing abundance, 





reaction cross-section, and average flux values, the equilibrium activity 
due to this reaction will be roughly 10% of the "el activity. Since in. - 
addition the delayed gamma from ofS is far less penetrating than that 
from ni6 it is reasonable to assume that the escaping ol3 delayed gamma 
activity will be negligdble. 
4. As a result of the above process of elimination, there remains only 
one induced activity in the coolant which need be considered; the ni6 
delayed gamma which results from the ol6 (n,p) ni6 reaction, The threshold 
for this reaction is 9.5 Mev. The reaction cross-section, averaged over 
the fission spectrum is 0.014x10-> barn. The half life of the nt product 
is 7.35 seconds, and the gamma energy is 6 Mev. Seventy five percent of 
the n6 delayed gamma decays are accompanied by this high energy gama. 
5. Based on the assumptions of articles 3 and 4, above, a rough and 
simple system was devised. The system consisted of 
(a) The Detector- An anthracene crystal and a 5819 
photomultiplier tube 
(b) The Indicator— An Electrometer, Model 4-565-1, 
Mfgr. Victoreen Instrument Company. 

The detector was heavily shielded to prevent any interference from 
stray radiations, and to give it directional properties. A check of the 
angularity at the detector, using a 77 mc Co®° source at a distance of 
one foot, indicated a beam width of approximately 9.54 degrees. 

The output of the detector was put to the Electrometer and null 
voltage readings used for plotting purposes. Later in the experimental 
work the detector output was fed to a Sanborn Recorder for continuous 


readings, and occasional null voltage readings were taken to correlate 





the Sanborn scale with the Electrometer readings for the previous runs. 
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Fig. 1 — Simplified Block Diagram of Experimental 


Recorder 


System Showing Relative Positions of Equipment 
With the physical layout of the system as shown above, and inasmuch 
as the only induced reaction of concern is the olé (n,p) ni6 reaction, 
the following constants may be assumed to hold for use in equation (19) 


of Chapter I: 


A = 1.70 Square Feet 

L = 26 Feet 

Q = 400 Cubic Feet 

T = £20 Seconds 

A = 2£§93 = 0,0943/Seconds 
Te SA) 





Equation (18) of Chapter I may now be written 
~ £9-0943N2Z0)(1.79)(26) 


R = KP e 400F 
_ 0- 2085 
or R = KP e = (3) 


In this form the equation will apply to the particular experimental 


system used. 


6. At this point, it is advisable to stop and consider the effect of 


flow on the equilibrium concentration of ni, 


From equation (10) of Chapter I 
No= N, agen *™) (4) 


8 oe 


where t = 0 at the effective time of entering the core. The number of 


activated nié 


nuclei in the coolant, at any time t while in the core, 
may be found from equation (4) above. 

Assuming that equilibrium is attained in passing through the core, 
the number of activated nuclei in the coolant, at any time t after leaving 
the core, may be obtained from the simple decay law (Equation (2) of 
Chapter I) as 

ee 2 ve (5) 

To determine the effect of flow on the equilibrium value of n26 

equations (4) and (5), above, must be PGnsitier ed simultaneously in the 


following manner: 


Productisoy Curve 
N« = Np. {i = ey 







Oecay Curve 
of + é t; Nz =ard en eC” 
7/772 

Let t) = loop time (or decay time before reentering the reactor) 


AG 


t, = Effective initial production time (or time which would have 
been needed to build up to same value of N, , along normal 
production curve, as that with which the coolant is reentering 


the reactor) 





ct 
is 


5 Effective total production time ™ ttt, 
t, = Actual time in core (actually a 4t) 
t. = Time from reactor outlet to detector 
At equilibrium; t) , tz , and t3 must be So related that 
No (lee > ) = Ng, (1-e7*"3) eA (6) 
Utilizing the fact that t) and t, can be expressed as "known" 
functions of flow and that t, = t,+t,, solving equation (6) for 


3 
e7 **2 one obtains 


reg: 
-aAt 7 
eure [Sam | (7) 
-@ 

Since 

eS er Alen tty) -dEe oae4 (8) 
then . anaes 

_at es -xr# 

e "3 E . as | e€ x (9) 

and Tee 2 tae | ue” ee ee (10) 
ae fs eg Alt +te 


Therefore at equilibrium, the number of activated nuclei per unit 


volume of coolant at the reactor outlet re 
¢ 


-rt 1-@e@ 
| No = No, (lee *)= Ng. [a J (11) 
This indicates that equation (3) of this chapter should be modified 
to include a factor for the effect of flow on the equilibriwm activity 


within the coolant. This equation then takes the form 
“rt, 


- l-e E 
2 Genesco ld (12) 
Substituting the constants from article 5, Chapter II into 


equation (12) 





ar - - 1 AGE HE (13) 
1.8 
c+ t. 2 = z ss Be | il Lae (14) 
C = 0.2085 
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and 1.63 
ae ee 
R = xe [ — . en (15) 





The total effect a flow on the activity indicated by the detection 
system, is best indicated by the following table. 


Table 1 — Correction Factor for Flow Effect 


Percent Varig bh ioe 


ie PP 22087 «Correction -*F Feeder Srecns 
1.00 0.947 0.812 0.769 == 
0.95 0.950 0.803 0.763 0.8 
0.90 0.953 0.793 0.755 1.8 
0.85 0.957 0,782 0.749 2.6 
0.80 0.961 0.770 0.740 3.8 
0.75 0.964 0.757 0.722 6.1 
0.70 0.968 0.742 0.718 6.6 
0.65 0.971 0.725 0.704 &. 5 
0.60 0.976 0.706 0.689 10.4 
0.55 0.981 0.685 0.672 2.1 
0.50 0.984 0.659 0.649 15.6 


Inasmuch as the dimensions and times used for calculations are 
only approximations, and estimates based on these approximations, the 
variations as a result of flow changes are not prohibitive. The actual 
dimensions and times involved in the STR power plant are such as to 


reduce this variation. 
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CHAPTER III 


EXPERIMENTAL RESULTS 


1, In Chapters I and II it has been shown that the induced activity in 
the reactor coolant is proportional to the power level of the reactor, 

and that it is possible to devise an instrument system which uses this 
activity as a means of power level determination. A system such as 
described in article 5, Chapter II, was installed at the Submarine Thermal 


Reactor Test Site and the following data obtained. 


Reading Equivalent 
% Power % Flow Reading "Scale" Corrected Reading 
for 


Background 
RUN NUMBER ONE 


12 


0.0 17.2 52.0 ——— 0.0 1.0.0 
5 17.2 60.3 —- Bias 8.3 
2.6 17,2 69.6 ——~ 17.6 17.6 
Bie 17.2 76.5 —_——~ 24.5 24.5 
9.5 17,2 85.0 ———— 33.0 33.0 
10.75 17.2 93.0 ae 41.0 41.0 
11.5 17.2 93.0 -_——- 41.0 41.0 
13.0 ae ary 97.8 ———— 45.8 45.8 
RUN NUMBER TWO 
2.9 96.5 101.0 ———— 53.0 56.5 
15.0 96.5 119.0 ——— 61.0 65.0 
16.0 96.5 122.0 — 64.0 68 .0 
18.0 96.5 129.0 ——— 71.0 1362 
21.5 96.5 143.0 —_—— 85.0 90.5 





Reading Equivalent 
% Power % Flow Reading "Scale" Corrected Reading 
for 


Background 
RUN NUMBER TWO (CONT. ) 


25.0 
27.5 
30.0 
33.0 
40.0 
41.0 
41.5 
42.0 
42.5 
44.0 
45.0 
46.0 
47.0 


47.5 


96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 


155.4 
155.0 
157.0 
172.7 
189.3 
192.8 
192.5 
197.0 
196.3 
202.0 
207.7 
209.0 
211.0 
Z15.0 
222.8 
224.4 
222.0 
226.0 
229.0 
225.0 
231.0 
231.0 


231.0 
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104.0 
103.0 
105.0 
122.0 
140.0 
144.0 
143.0 
148.0 
147.0 
153.0 
159.0 
161.0 
163.0 
168.0 
175.0 
177.0 
175.0 
179.0 
182.0 
178.0 
184.0 
184.0 


184.0 


° ~~ 





Reading Equivalent 
% Power % Flow Reading "Scale" Corrected Reading 
for 


Background 
RUN NUMBER TWO (CONT.) 


55.0 


96.5 


96.5 


235.0 
245.0 
247.0 
253.0 
255.0 
259.0 
261.0 
263.0 
265.0 
268.5 
267.0 
279.0 
285.0 
291.0 
277.0 
293.5 
294.3 
292.8 
294.0 
301.0 
308.0 
308 .0 


310.0 
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177.0 
187.0 
189.0 
195.0 
197.0 
201.0 
203.0 
205.0 
207.0 
210.5 
209.0 
221.0 
227.0 
233.0 
219.0 
235.5 
236.3 
234,8 
236.0 
243.0 
250.0 
250.0 


252.0 


188.0 
199.0 
201.0 
208 .0 
210.0 
214.0 
216.0 
218.0 
220.0 
224.0 
222.0 
235.0 
242.0 
248 .0 
233.0 
245.0 
252.0 
250.0 
251.0 
258.0 
266.0 
266.0 


268 .0 





Reading Equivalent 
% Power % Flow Reading "Scale" Corrected Reading 
for 


Background 
RUN NUMBER THREE 


37.5 
38.0 
40.0 
aie 
67.5 


69.0 


5<0 
13,5 
25.0 
58.0 
oe eS 
59.0 


61.0 


5.0 
10.0 
15.0 
20.0 
25.0 
30.0 
35.0 


40.0 


80.0 
80.0 
80.0 
80.0 
80.0 


80.0 


RUN NUMBER FOUR 


96.5 
96.5 
96.5 
96.5 
96.5 
96.5 


96.5 


RUN NUMBER FIVE 


96.5 
96,5 
96.5 
96.5 
96.5 
96.5 
96.5 


96.5 


187.0 
207.0 
207.0 
227.0 
295.0 


301.0 


51.0 
84.0 
113.0 
212.0 
212.0 
210.0 


237.0 


57.0 
73.0 
88 .0 
101.0 
115.0 
124.0 
139.0 


143.0 


15 


135.0 
155.0 
155.0 
175.0 
243.0 


249.0 


17.0 
50.0 
79.0 
178.0 
178.0 
176.0 


203.0 


13.0 
29.0 
44.0 
67.0 
71.0 
80.0 
95.0 


Sele 





Reading Equivalent 
% Power % Flow Reading "Scale" Corrected Reading 
for 


Background 
RUN NUMBER FIVE (CONT.) 


45.0 
50.0 
55.0 


60.0 


RUN NUMBER SIX 


8.0 
10.0 
10.5 
15.0 
20.0 
25.0 
29.0 
30.0 
34.0 
35.0 
40.0 
44.0 
45.0 
48.0 


49.5 


96.5 
96.5 
96.5 


96.5 


92.5 
96.5 
96.5 
96.5 
9G, 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 


96.5 


156.0 
170.0 
183.0 


201.0 


770.0 
74.0 
i330 
81,0 

100.0 

110.0 

122.0 

125.0 

138.0 

138.0 

144.0 

161,0 

169.0 

177.0 


174,0 


16 


1140 
1250 
1350 


1480 


490 
520 
510 
. 570 
710 
790 
890 
900 
990 
990 
1040 
1170 
1220 
1270 


1200 


122.0 
126.0 
139.0 


157.0 


20.0 
24.0 
23.0 
31.0 
50.0 
60.0 
72.0 
75.0 
88.0 
88.0 
94.0 
111.0 
119.0 
127.0 


124.0 


152.0 
170.0 
188.0 


212.0 


28.0 
34.0 
32.5 
44.0 
70.5 
84.5 

102.0 

106.0 

124.0 

124.0 

133.0 

157.0 

168.0 

179.0 


175.0 





Reading Equivalent 
% Power % Flow Reading "Scale! Corrected Reading 
for 


Background 
RUN NUMBER SEVEN 


5.5 
10.0 
15.0 
20.0 
25.0 
30.0 
34.5 
35.0 
39.5 
40.0 
41.0 
44.0 
45.0 
49.0 
49.5 
50.0 
50.5 
55.0 
60.0 
65.0 
67.5 
70.0 


75.0 


96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 
96.5 


96.5 


52.0 

60.0 

75.0 

86.0 
104.0 
118.0 
127.0 
127.0 
144.0 
144.0 
152.0 
155.0 
155.0 
168.0 
163.0 
162.0 
175.0 
166.0 
173.0 
209.0 
218.0 
223.0 


236.0 


17 


310 
350 
460 
530 
650 
740 
790 
800 
910 
910 
970 
‘980 
980 
1060 
1030 
1030 
1100 
1050 
1100 
1320 
1400 
1420 
1500 


14.0 
22.0 
37.0 
48.0 
66.0 
80.0 
85.0 
85.0 
106.0 
106.0 
114.0 
117.0 
117.0 
130.0 
125.0 
124.0 
137.0 
128.0 
135.0 
171.0 
180.0 
185.0 


198.0 





% Power 


% Flow 


Reading 


"Scale" 


Reading 
Corrected 


.for 


Equivalent 
Reading 


Background 
RUN NUMBER EIGHT 


40.0 
45.0 
50.0 
55.0 
60.0 


65.0 


96.5 
96.5 
96.5 
96.5 
96.5 


96.5 


RUN NUMBER NINE 


2. 


70.0 
75.0 
80.0 
85.0 
90.0 
95.0 


100.0 


96.5 
96.5 
96.5 
96.5 
96.0 
96.5 


96.5 


145.0 
163.0 
170.0 
176.0 
184.0 


191.0 


203.0 
220.0 
233.0 
246.0 
272.0 
290.0 


301.0 


860 
970 
1010 
1050 
1100 
1140 


1220 
1320 
1400 
1480 
1640 
1735 
1815 


47.5 
65.5 
72.5 
78.5 
86.5 


9355 


\ 
228 .0 


245.0 
258 .0 
271.0 
297.0 
315.0 


326.0 


116.0 
160.0 
177.0 
192.0 
212.0 


228.0 


249.0 
266.0 
280.0 
294.0 
322.0 
342.0 


354.0 


A plot of the actual reading versus percent power (Figures fwo and 


three) indicates that the induced activity in the coolant is, indeed, 


linear with respect to power. 


However, there is an obvious variation in 


the slope of the line for each run; and a different background constant 


for each run, 


equations for the various lines are: 


AsAsoceach line can be represented by the equatiosm R,= m, P+ C, the 
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Run 


Number 


One}; 


354 8 + 58 


1 
Run Number Two; R, = 333 P + 52 
Run Number Three}; R3 360P + 52 
Run Number Fours Rg © 332P + 34 
Run Number Fives Re = 262P + 44 
Run Number Six; Re 251P + 50 
Run Number Seven; Ry 265P + 38 
Run Number Eights; Rg w= 145P + 97.5 
Run Number Nine} Ro 326, 5P-25 


In order 


for background, and then corrected for the differences in slopes. 


to compare the various runs, each run must be corrected 


Con- 


verting all data points to equivalent readings, which will correspond to 


run number one, 


involves the following mathematical processes with each 


data point: 
RY = Reading Corrected for Background = Rn - CG, 
Ry, — Equivalent Reading - a R. 


The results of this manipulation are tabulated in the data tables, 


and are plotted in Figure Four. To briefly describe Figure Four, the 
resulting relationship between the equivalent reading and percent power 
is linear. Of the 135 experimental points, 102 (or 75.6% of the total) 
vary by four percent or less from the value which would be predicted from 
this line. Of this 102, it is interesting to note, that 40 (or 29.6% of 
the total) vary by approximately zero percent. It should also be noted 
that the data has not been corrected for flow effect, thus indicating that 


the effect of flow is not as great as that determined theoretically in 


r 


Chapter II. 
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It may be reasoned that the variations in the slopes of the plots 
of original readings are due to impurities and flow rather than any 
effect of control rod configuration as; 

(a) there was little change in the control rod configuration 
during the caper edited periods and 

(bo) correcting the experimental readings for batkground, 
plotting these, and extrapolating all runs to zero power 
will result in all runs intersecting at a common point at 
zero power. 

Theoretically any radial effect of control rod configuration will be 
well masked and eliminated by mixing in the header prior to entering the 
external loop piping. Changes in control rod configuration may have some 
influence on the coolant activity at the reactor outlet due to axial 
displacement of the point of maximum flux within the core. However, no 
calculations can be made on this effect without going into the actual 
Biristruction and loading of the core. Since, for normal operation, con- 
trol rod configuration should be such as to give as uniform a power and 
flux distribution as possible, the flux peak will not be too pronounced 
and will not drift much during the normal operating period between thermal 
calibrations, This effect may be such as to either increase or decrease 
the apparent power level depending on the reactor design. It is estimated 
that, over the normal operating life of the core, this effect will be of 
such smal] magnitude as to be easily accounted for during periodic thermal 


calibrations of the entire instrument systen. 
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CHAPTER FOUR 


EVALUATION OF SYSTEM 


1. It ms been shown that the number of activated nuclei per unit volume 
of coolant, at the primary coolant outlet of a given reactor system, can 
be described by the sae sag 

N oo 7 Paes (= | 

or No? <= (sa | (1) 


From the simple radioactive decay law, the decay rate of this unit 


volume at any time t. after leaving the reactor is 





fe en + - re 
tt = -ANg = -£¢.. [ ; atta | e J (2) 
As flux is directly proportional to power, this equation can be 
written -ré 
dNe K Sy Pa C /-e < ] aes 
Jt a Sie.) Fa ie er NE, 4F 4) (3) 


A directional gamma detector may be installed at some point on the 
external coolant piping, and an indication obtained which is proportional 
to the power level of the reactor, This indication would have the advan- 
tages, over conventional neutron detection systems, of being independent 
of the control rod configuration within the reactor core; and would use 
crystalline materials in which there is no apparent permanent damage due 
to radiation (i.e., no aging). 

2. While the experimental data does not support this conclusion as fully 
as might be desired, it is actually quite encouraging considering that: 
(a) Data was taken on an tas available, not to interfere’ 
basis, and data on characteristics which it would have 


been desirable to investigate was not available; 
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(b) The general condition at the reactor system, as a 
result of previous and current tests, was far below 
standard; and often the operating wmnditions imposed 
were abnormal 3 
(c) The educational background of the writer, at the time 
the experimental work was done, was imperfect (to say 
the least) in this field. As a result, it was not 
realized that many other points should be considered 
(e.g., concentration of impurities in coolant, etc.). 
3. The dependence of the indication upon the flow rate of the coolant 
for an existing reactor system can be minimized quite easily by Borer 
physical location of the detector. As an example; the percent setateuea 
in the flow correction factor for the experimental system was 15.6 over 
a flow range of 50-100 percent. Had the distance from the reactor outlet 
to the detector been only one fourth as much (i.e. 6.5 feet instead of 
26 feet) the correction factors would have been 
for F => 1.003 Correction factor = 0.947e  =0.898, 
and for F = 0.50; Correction factor = 0.984e  =0.886; 
or a percent variation in the correction factor of 1.3 over a flow range 
of 50-100 percent. 
For an initial design, the flow correction factor in equation (3) 
- 0.69 5a = 


can be written as [ + -™ "7 -0.693a,5 
1 7 693 7, © la 


Where T = Time for a complete flow cycle; 
1 Half life of the radioactive species of concern}; 
“,= Fraction of T required to traverse the core} 


“= Fraction of T required to reach the detector. 
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Examination of the factor in this form indicates that for minimum effect 
of flow on the indication: 
(a) The ratio of flow cycle time to half life should be 
as large as possible; 
(b) The fraction of flow cycle time spent in the core 
should approach one, as closely as possible and 
(c) The fraction of flow cycle time to reach the detector 
should approach zero, as closely as possible. 

The major effect of (a) and (b) will be to influence the equilibrium 
concentration of the activated species in the coolant. From the design 
point of view, (c) is the most critical; yet the simplest to control. 

Based on the design form of the correction factor, it is possible 
to plot design curves of percent variation in correction factor (effec- 
tively percent error in indication) versus the ratio of flow cycle time 
to half life. A different curve will be obtained for each combination 


of «, and & however, the only combinations which appear to be of 


2? 
any practical interest are those of 0.54%,£1.0 and 0.04%,¢0.5. 
Plotting these design curves, which involves a multitude of simple 2 
calculations, is considered by the writer to be beyond the scope of 

this paper. Through the use of the design curves or the design form 

of the correction factor, it is possible to optimize the position of 

the detector and the ratio of loop time to core time for any given system, 
using any desired coolant, so that the maximum percent error in the in- 
dication will be within any desired range. 

4. There are two effects on the gamma activity in the reactor coolant 


which are reflections of the chemical elements which make up the coolant 


and the impurities present. 
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The first effect is on the decay constant, A . The A used in the 
equation which describes the induced delayed gamma activity in the reactor 
coolant must be a weighted average A. As such, it will be relatively 
strom ly influenced by the presence of small amounts of other than normal 
radioactive substances with A's which vary considerably from the average 
for the pure coolant. 

The second effect is on the slope of the relationship between the 
power level and the gamma indication. The slope of this relationship is 
directly proportional to the macroscopic cross-section of the reactor 
coolant. Obviously this Chitediamat ones) be strongly influenced by 
small amounts of large cross~—section elements, either normally present 
in the coolant or introduced as impurities. 

If the indication of power level is based on the gross gamma activity, 
as was the case with the experimental system, slight changes in the con- 
centration of impurities present, or the introduction of very small amounts 
of fission products into the coolant system, can produce marked changes in 
the slope of the relationship and the effect of flow on the indication. 
While no experimental data can be produced to support this, as data was 
not available on specific coolant composition during the tests, it was 
known that the composition was variable from operating period to operating 
period. This would account for the variation of slope from run to run in 
the experimental data, and for some of the seemingly "bad" points on the 
resultant plot. 

The first of these effects can be effectively eliminated by a rather 
simple device. If only one specific activity were to be observed, and all 
other gamma activity eliminated from the detection system, the A of con- 


cern would be the A for the specific activity and would remain constant, 
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Therefore a simple fixed range pulse height selector placed in the circuit 
ahead of the indicator would pass only the scintillations produced by the 
gammas from the selected decay process, and the power level indication 
would depend only on the selected activity. This also suggests the pos- 
sibility of introducing small amounts of some tracer element into a 
coolant which did not have a neutron induced delayed gamma reaction of 
its own which would be suitable for detection purposes. 

The second effect is more difficult to control, however it should 
be remembered that under normal operating conditions the concentration 
and number of impurities in a reactor system would be much smaller than 
in the experimental or prototype system. The control of impurities in 
water and. other liquids is not a new problem to engineers and can often 
be handled in a relatively simple manner. In a pressurized water system , 
this could be accomplished by partial continuous or intermittent circu- 
lation of the coolant through a mechanical filter, then resin demineralizers. 
It may be that both cation and anion exchange should be utilized for maxi- 
mum effectiveness. This would be a system similar to some of those proposed 
for use in connection with boiling-water reactor systems. As the system 
would be a closed system, the impurities would not built up rapidly, and 
continuous checking of the impurity concentration would not be required. 
Intermittent checks of the impurity concentration could be accomplished by 
a modification of the low level monitor system developed by E. I. duPont 
de Nemours and Company (2). Chemical analysis could of course be used, 
but is somewhat more involved, and requires personnel trained in chemical 
analysis techniques. With controlled impurity conventration, no slope 


correction circuit would be required in the metering circuit, but could 
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be included as an added precaution in the event that an occasional thermal 
calibration should indicate a slope correction was needed. 
5. It is the opinion of the writer that it is quite possible to design 
a system for power level determination which uses the neutron induced 
gamma activity in a reactor coolant as the source for the indication. 
Such a system can be made as accurate as desired, and independent of such 
things as flow rates, control rod configurations, deterioration of materials 
(aging), and impurities in the coolant. By suitable choice of certain 
design variables, the system can be made applicable to any existing or 
proposed nuclear reactor system which circulates a coolant through a 
region of neutron flux. The system may be designed for any type of 
coolant, organic or in-organics and is equally applicable to boiling 
reactor systems. In the latter, it appears that the impurity problem is 
not as serious as in those systems where no phase change (i.e. liquid- 
vapor) in the coolant is encountered. 

No attempt has been made to evaluate the transient behavior of this 
system, however observations of the experimental system during the experi- 
mental period indicated that it does have many desirable features under 


transient conditions. 
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